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EFFECTS  OF  ION-ACOUSTIC  INSTABILITY  ON  LIGHT 
ION  BEAM  TRANSPORT  IN  DEUTERIUM  CHANNELS 


I.  Introduction 

In  a  previous  work^-,  the  aagnetohydrodynamic  response  of  plasma  channels 
to  propagating  Ion  beams  has  been  Investigated.  It  was  found  that  3  MeV,  0.5 
MA,  50  nsec  proton  beams  could  be  transported  In  one  cm  deuterium  plasma 
channels  without  major  disruptions  of  the  channel.  This  conclusion  was 
further  supported  by  the  results  of  a  review  on  instablitles  of  ion  beam 
propagating  in  preformed  channels  .  However,  when  the  primary  ion  beam 
current  is  increased  to  several  MA,  Fig.  3  in  ref.  1  for  example  shows  that 
large  values  of  Te/T^  can  be  reached  in  the  channel  (of  order  10)  favoring  the 
growth  of  ion-acoustic  instabilities  and  that  unfavorable  structure  in  the 
radial  B  field  profile  may  appear,  raising  questions  about  the  confinement  of 
injected  ions  in  the  channel.  In  both  of  these  cases,  classical  transport 
which  was  assumed  previously  is  not  a  good  assumption. 

In  this  work,  we  modify  the  transport  coefficients  due  to  the  onset  of 
ion  acoustic  instabilities,  we  determine  the  conditions  under  which  ion 
acoustic  instabilities  are  most  likely  to  occur  and  recalculate  the  channel 
response  when  the  effects  of  these  instabilities  are  included.  Changes  in 
channel  resistivity  modify  axial  energy  losses  of  the  beam  during  propagation 
and  consequently  the  optimum  channel  density  for  minimizing  these  energy 
losses  .  In  a  first  section,  we  look  at  the  ion-acoustic  instability  and  its 
consequences  on  the  physics  of  the  problem.  Then  in  a  second  section,  we 
present  the  results  of  this  study  and  we  draw  conclusions.  It  appears  as 
though  ion-acoustic  instability  and  the  associated  anomalous  beam  slowing  down 
are  not  very  important  effects  on  the  ion  beam  propagation  in  the  plasma 
channel. 
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II.  Ion-acoustic  Instability 


Ue  have  already  described  the  treatment  of  ion-acoustic  instabilities  in 
fluid  calculations^.  For  sake  of  completion,  we  outline  this  treatment  here, 
but  the  emphasis  will  be  on  the  differences  between  the  present  case  and  the 
previous  one.  In  the  present  case,  a  return  current  flows  through  the 
channel,  heating  up  the  electrons.  If  we  denote  the  current  drift  velocity  by 
VD,  the  lon-acoustlc  instability  is  excited  when 
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where  V  is  the  ion  sound  speed  given  by 
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All  the  symbols  are  standard  (M  is  the  ion  mass,  m  is  the  electron  mass  for 
instance)  and  the  temperature  is  expressed  in  ergs.  The  expression  for  the 
instability  threshold  is  valid  when  ZTe  »  T^.  When  the  Debye  length  XDe  is 
substituted  by  its  expression  as  a  function  of  density  and  electron 
temperature,  Eq.  (1)  can  be  expressed  as  a  condition  for  VD/Ve  where  Ve  is  the 
electron  thermal  velocity.  For  a  deuterium  plasma  (Z  •  1,  M/m  •  3700),  VD/Ve 
has  a  maximum  around  T^/Tg  ■  0.45.  For  larger  values  of  T^/Tg  -  this 
situation  occurs  in  the  channel  at  its  edge  because  of  shock  heating  -  we  use 
the  following  expression  which  is  the  approximate  condition  for  the  two-stream 
instability 


V_/V  -0.84  ==-+0.2. 
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A  smooth  transition  between  the  regimes  defined  by  relation  (1)  and  (3)  is 
assumed  and  has  been  Implemented  into  the  code. 

Once  the  threshold  conditions  have  been  met  the  effects  of  lon-acoustlc 
instabilities  on  the  plasma  behavior  take  place.  These  effects  have  been 
discussed  in  ref.  4  and  arise  from  an  anomalous  collision  frequency,  anomalous 
electron-ion  energy  exchange  and  Increased  electric  field  which  generates  the 
return  current.  Of  particular  importance  to  the  present  work  is  the  enhanced 
electric  field  which  decelerates  the  beam  ions  in  the  channel  and  thus 
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Increases  their  slowing  down.  This  electric  field  Is  simply  calculated  as 
where  the  anomalous  collision  frequency  has  been  chosen 


v  -  10"3  u>  . 
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A  factor  of  10  Instead  of  10  was  also  tried  In  the  computation,  but  the 
difference  was  not  significant.  Anomalous  heat  transport  was  not  Included. 


III.  Results 


In  this  section,  we  first  look  et  the  physical  regime  for  which  the  ion- 
acoustic  instability  is  important.  From  the  results  of  ref.  1.  we  have 
already  seen  that  the  the  temperature  ratio  T^/Te  in  the  channel  at  the  end  of 
the  beam  pulse  varies  between  0.7  for  a  400  kA  proton  beam  and  0.1  -  0.2  for 
the  2  MA  case.  From  the  threshold  condition,  we  know  that  the  instability 
will  occur  when  VD/Ve  lies  around  0.8  -  0.4  for  the  temperature  ratios  just 
quoted.  Taking  the  larger  value,  expressing  Ve  as  a  function  of  Te  and  using 
the  relation  between  the  temperature  and  current  density  found  in  ref.  1,  we 
find  that  the  equivalent  threshold  condition  which  must  be  satisfied  by  the 
current  density  for  the  instability  to  take  place  is 
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where  j^  is  in  A/cm  ,  ne  in  cm  and  e  the  ion  beam  energy  in  MeV.  This 
expression  is  valid  for  a  proton  beam  in  a  deuterium  background  channel  for  a 
50  nsec  ion  beam  pulse.  In  the  bunched  stage,  we  know  that  jj,  t  is  constant 
and  condition  (5)  becomes  more  easily  satisfied,  especially  when  the  numerical 

IQ 

factor  on  the  right  hand  side  drops  to  1  for  a  10  nsec  pulse.  For  ne  «  10iO 
cm”^,  e  •  1 


Jb  >  107  A/ cm2 

in  the  bunched  stage.  Another  factor  which  helps  in  meeting  that  criterion  is 
the  reduced  channel  density  after  the  beam  has  been  Injected.  This  reduction 
factor  has  been  estimated  in  ref.  3  and  should  also  be  taken  into  account. 

We  now  present  results  for  a  bunched  case,  j^  ■  107  A/cm^,  hq  *  4  10*^ 
cm  and  e  ■  2  MeV.  Profiles  for  density  and  temperature  as  a  function  of 
channel  radius  with  and  without  the  instability  are  shown  near  the  end  of  the 
pulse  at  t  "  8.2  nsec  in  Figs.  1  and  2.  For  the  anomalous  transport  case, 
heating  is  higher  at  the  channel-gas  blanket  Interface.  Also,  because  of  the 
larger  current  resulting  from  the  anomalous  resistivity,  the  B-field  has  more 
structure  than  in  the  classical  case. 

As  mentioned  before,  one  of  the  most  important  quantities  which  comes  out 
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of  this  study  Is  the  decelerating  axial  electric  field  in  the  channel.  The 
slowlng-down  of  the  injected  beam  is  expected  to  be  larger  in  the  unstable 
case  because  of  the  anomalous  resistivity.  The  electric  field  shown  in  Fig.  3 
has  been  averaged  over  radius  according  to  the  expression 

fE  (r)  2xr  dr 
<  E  >  -  —2 - 2 - 

irr  , 
ch 

and  includes  also  the  colllsional  slowing-down  of  the  beam  in  the  channel.  It 
is  shown  as  a  function  of  initial  channel  density  for  the  same  conditions  as 
indicated  before  at  t  ■  8.2  nsec,  i.e.,  just  before  the  end  of  the  beam 
pulse.  From  this  figure,  we  see  that  the  average  electric  field  for  the 
unstable  case  is  larger  than  for  the  classical  one  but  not  by  a  large 
amount.  The  difference  goes  from  about  one  percent  at  n^  »  10  °  cm  to  45Z 
at  n^  *  2.5  10  cm  .  The  field  reaches  a  minimum  as  found  in  ref.  3  when 
energy  losses  due  to  motion  of  the  channel  (V  x  B  term  in  the  electric  field 
inversely  proportional  to  density)  equal  losses  due  to  beam  deposition 
(proportional  to  density).  The  density  value  at  which  this  minimum  occurs  can 
be  compared  with  the  simple  estimates  given  in  ref.  3.  For  the  present 
parameters,  these  estimates  give  the  minimum  in  energy  losses  for  n^  -  1.8 
101®  cm“^  whereas  the  numerical  results  in  Fig.  3  show  it  at  101®  cm-"*.  For 
the  most  efficient  transport,  channels  should  operate  on  the  high-density  side 
of  the  minimum  (since  the  losses  do  not  go  up  as  fast  with  a  variation  in 
density  than  on  the  other  side).  It  is  also  the  side  where  effects  of  ion- 
acoustic  instability  are  smaller.  For  this  reason,  it  does  not  seem  that 
anomalous  transport  plays  an  important  role  in  beam  transport  studies. 

In  summary,  we  have  looked  at  the  effects  of  ion-acoustic  instability  on 
channel  transport.  Anomalous  transport  coefficients  have  been  Included  in  a 
one-dimensional  radial  mhd  code.  The  effects  of  anomalous  transport  have  been 
shown  to  be  significant  only  in  the  bunched  stage  when  the  initial  channel 
density  is  below  its  optimum  value.  Because  of  these  characteristics,  it  is 
our  conclusion  that  effects  of  ion-acoustic  instability  do  not  seem  to  be 
important  for  beam  transport  studies. 


Plasma  channel  density  profiles  with  and  without  ion-acoustic 
instability  for  »  107  A/cm3,  nQ  ■  4  x  1037  cm  3  and 
e.  «  2  MeV  at  t  »  8.2  nsec. 
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Fig.  2. 


Plasma  channel  temperature  profiles  with  and  without  ion-acoustic 
instability  for  same  parameters  as  Fig.  1. 
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Fig.  3. 
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Equivalent  axial  electric  field  with  and  without  ion-acoustic 
instability  as  a  function  of  initial  channel  density  at  t  ■  8  nsec 
for  •  10  7  A/co2  and  ■  2  MeV. 
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